Physical realization for Riemann zeros from black hole physics 
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According to a conjecture attributed to Polya and Hilbert, there is a self-adjoint operator whose 
eigenvalues are the the nontrivial zeros of the Riemann zeta function. We show that the near-horizon 
dynamics of a massive scalar field in the Schwarzscild black hole spacetime, under a reasonable 
boundary condition, gives rise to normal mode frequencies that coincide with the nontrivial Riemann 
zeros. In achieving this result, we exploit the Bekenstein conjecture of black hole area quantization, 
and argue that it is responsible for the breaking of the continuous scale symmetry of the near horizon 
dynamics into a discrete one. 
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/. Introduction. — The Riemann zeta function is of 
great importance in number theory, and is of interest to 
physicists because of its connection to quantum chaos 
and random matrix theory [H-Q- Starting with the 
Dirichlet series 



oo 



(1) 



n=l 



vifhich converges for TZe s > 1, Riemann [4] analytically 
continued it over the whole complex plane by defining an 
appropriate contour integral that reduced to Eq.(IT]) for 
TZe s > 1. Writing s = (si + is2), the only pole is on 
the real axis at si = 1, and the "trivial zeros" are on the 
real axis at negative even integer values of si. Accord- 
ing to Riemann hypothesis, the only other zeros of C(s)i 
infinite in number, lie in complex conjugate pairs on the 
line si — 1/2 in the complex s-plane. This has been nu- 
merically tested over extensive samples by Odlyzko Q, 
but has not been proven mathematically [7]. According 
to a conjecture attributed to Polya and Hilbert, these ze- 
ros on the upper half of the complex plane may be the 
real eigenvalues of a self-adjoint operator in quantum me- 
chanics. Previous work on this topic are due to, amongst 
others, Hua Wu and Sprung Berry and Keating 2], 
Connes [9i] , and Sierra and Townsend [13] ■ 

The purpose of this paper is to suggest such an example 
from black-hole (BH) physics. Specifically, we consider 
the Klein-Gordon spectrum of a massive scalar field just 
outside the horizon of a Schwarzschild BH [ll| . It is well 
known that the near-horizon dynamics of the scalar field 
is described by a scale-invariant inverse-square potential 
with a coupling strength that is a function of the eigenval- 
ues of the scalar field [l^ . This is quite different from a 
nonrelativistic self-adjoint operator with a constant cou- 
pling. The scale invariance is broken by other terms in 
the Hamiltonian at larger distances from the horizon. 
The correspondance of the eigenvalues of the scalar field 



with the Riemann zeros depends crucially on the assump- 
tion of scale invariance with energy-dependent coupling, 
and therefore refers to only a subset of the deeply bound 
eigenvalues that come from the near-horizon dynamics. 
There are two other important ingredients in obtaining 
the desired result. The inverse square potential 's cut 
off by a brick-wall as introduced by 't Hooft [ll| very 
near the horizon, but we demand that the final spec- 
trum should not depend explicitly on the brick-wall ra- 
dius. Lastly, we make the crucial assumption that the 
area of the BH is quantized in units of the Planck length 
squared. In order to preserve this quantisation under 
a scale transformation in D-spatial dimensions {D > 1), 
the continuous scale symmetry breaks into a discrete one. 
In section II we analyze the massive scalar field equation 
in Schwarzschild black hole spacetime. In section HI the 
near-horizon approximation is made that gives a scale 
invariant Hamiltonian. The normal modes which have 
been found are identified with the nontrivial Riemann 
zeros of the zeta function. Finally we conclude in section 
IV with a discussion. 

//. Scalar field in Schwarzschild black hole spaceime. — 
As mentioned in the introduction, we go beyond the non- 
relativistic self-adjoint operator to obtain our desired re- 
sult. We begin our discussion with Schwarzschild black 
hole, which is a (Z? -I- l)-dimensional spacetime obtained 
from the vacuum Einstein Equation. For simplicity we 
consider a (3 + l)-dimensional spacetime defined by the 
metric 



ds^ 



-Ndt^ + N-^dr^ + r^dW 



where 
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M is the mass of the black hole and dfl^ is the metric 
over 5^. The dynamics of a massive scalar field of mass 
TOs in this black hole spacetime is given by the Klein 
Gordon equation 
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Unlike the BTZ black hole [Tl,[i3, this field equation is 
hard to solve exactly. However in order to get an idea 
of the normal modes for the scalar field, a careful look is 
needed for the corresponding near-horizon approximation 
of the eigenvalue equation. Because of spherical symme- 
try, we can write the ^-th partial wave-function for the 
normal mode w as 



1 



-iuJi t 



Yin 



». (5) 



In the above, r+ = 1M is the radius of the horizon, and 
a;; the eigenenergies of the scalar field. Substituting the 
ansatz in Q the radial eigenvalue equation is obtained 
to be 



1 



BH- 



The operator ilsH hi^ explicit form is written as 



(6) 



BH 



r^N — 
dr 
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dr 
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Henceforth we consider only the I = eigenstates that 
are non-degenerate, and arc relevant for our purpose. For 
simplicity, the / = quantum number is supressed from 
now on. Since ^ cannot be solved exactly, we proceed to 
obtain some information from the near-horizon approxi- 
mation. In the spirit of 't Hooft 11], to avoid divergences 
in free energy at the horizon, we apply a "brick wall" 
boundary condition that forces the scalar field to go to 
zero at the radius r = r+ + th where r+ + rn is slightly 
larger than the horizon radius r+ = 2M . Our final result 
is independent of rn- 

III. Near-horizon dynamics. — In this section we 
consider near-horizon dynamics of a scalar field in 
Schwarzschild black hole spacetime. In terms of the near- 
horizon coordinate z = (r — r_|_), the radial equation in 
the near-horizon limit is given by 



Hbh^uj{z) = 0, 



where 



H 



BH 



dz^ 



(8) 
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Terms of 0{l/z) or lower have been neglected in this 
limit. In general the operator like © occurs in many 
branches of non-relativistic physics also, but the crucial 
difference is the presence of the scalar field frequency lo 
as the coefficient of inverse square term. This will in fact 
help us achieve our goal. In order to show the relevance of 
to the Riemann zeros, let us now calculate quantum 
mechanical spectrum of ([8]). There are two independent 
short-distance solutions of the form 



(10) 



where 



and 



E = 2Muj 



(11) 



(12) 



The most general short-distance solution is a linear sum 
of both, 



(13) 



where Ai and A2 are constants. Note, however, that 
Hbh given by Eq.® is a scale invariant operator and in 
general it does not have any bound state solutions in the 
interval z G [0, 00] because the solutions are not normal- 
izable on the half line. However the asymptotic behavior 
of the radial equation ^ will make Fuj{z) normalizable 
at large distance. Restricting the dynamics in a finite 
interval may therefore give rise to bound state solutions. 
In the interval z S [rn , L] the operator ([9]) will be sym- 
metric in a domain D(Hbh) if it satisfies the condition 

^{z)*HBH(l>{z)dz^ ( (HBHiAz))* 4>{z)dz 
= ^{Ly<i>'{L)-^'{Lyc^{L) 

- i^{rH)*4>'{rH)+^'{rH)*^{rH)^Q, (14) 

where ■0(2): 0(^) G D{Hbh)- We have assumed that the 
field goes to zero asymptotically at a large distance ~ 
L. It can be easily seen in the following way. For large 
distance the effective radial differential equation, after 
neglecting terms of 0{1/ z) and lower in Eq.Q, becomes 



dz^ 



{ml - uj'^)F^{z) , 



whose decaying solution is of the form 

T^{z) = e-^- 



(15) 



(16) 



Note that the quantity (m| — uP') under square root 
is positive for bound state solutions, and L"^ ~ 
\J (m| — up-'). Therefore the condition (|14p reduces to 

^p{z)*HBH(f>{z)dz - / (Hbhi/Kz)) (j>{z)dz 

TH Jth 

= V'(^ff)Xrff)-V^(rff)*0'(rH) = 0, (17) 

which depends only on short distance cutoff r^j- The 
boundary condition consistent with (fT7|) will then be 



D[HBH)^{Hz)-4'{rH)^0}, 



(18) 



which makes the operator self-adjoint too |15|]. As men- 
tioned above, the field action corresponding to the op- 
erator ^ is scale invariant under the transformation 
t — > a^t and z — az, where a is the scale factor. The 
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scale invariance implies that if x(z) is an eigen-state of 
the operator Hbh with eigenvalue cbh then Xa = x(f ) 
is also an eigen-state of the same operator with eigen- 
value ^fir-- Note that the eigenvalue esH has nothing to 
do with the scalar field eigenvalue ui. In general a is a 
continuous parameter, which means the negative energy 
ground state will become unstable for a — > 0. One way to 
avoid this problem is to ask for the possible self-adjoint 
extensions which then can make the ground state 
stable. The self-adjoint extension parameter is related 
to the boundary condition. In our problem of black hole 
spacetime, however, black hole area quantization from 
quantum gravity effects may eventually break the contin- 
uous scale symmetry of the problem. As conjectured by 
Bekenstein [3, the black hole horizon area A — \i:r\ 
has a discrete spectrum 

= e^n-, n e N+ , (19) 

when it is treated as an operator in quantum gravity 
theory. In Eg. p^ above, we have used units G = c—\, 
so that the Planck length Ip ~ fi,^/^, and the Planck area 
is proportional to ft. The quantity eft is the constant 
spacing between the consecutive levels. The radius of the 
horizon goes like tj^^^ ^Jn through area quantization for 
D = 3. We exploit this quantization while making the 
scale transformation 

az — f ~ r+„j — ar — ar+„ . (20) 

From the second and third equality we see multiplying a 
to r+jj will take to a different only if a is quantized, 
implying that scale symmetry is broken into discrete scale 
symmetry by the area quantization. In general (|20p gives 
rise to the condition 

Q!m = v^,"ieN+ (21) 

The above scale transformation is rather ugly because 
it is different for each spatial dimension D. It is more 
satisfactory, instead, to take the scale transformation to 
be 

am=m,meN+ (22) 

which is valid for general space dimension D. 

Let us now consider the exact differential equation for 
the scalar field given by ([6]). Solutions of this equation 
which are globally valid are not scale invariant. How- 
ever the near horizon limits of the solutions are certainly 
scale invariant. Suppose J-\u)(^z) and J^2tj(z) are two in- 
dependent solutions of the differential equation (jS]). The 
general solution would be of the form 

T^{z) = BiFi^[z) + B2F2^{z) . (23) 

Due to short distance scale invariance we can construct 
a set of wavefunctions of the form Ti^{a^z)^am G 



Although Jlj(a~-^z) are not solutions of ([6]) we know that 
the short distance limit IiuIz^q ^^^{a^ z) ~ Ti^{a~^z) 
are solutions near the event horizon. Note that usually 
the scale transformation z ^ az takes the wave-function 
function il){z) — ilj{z/a). We then take these solutions 
and make a linear sum of them to impose a symmetric 
boundary condition of the form (|18p. which implies 

oo 

^ J-.(a„iz) = 0, (24) 

a™ = l 

at the brick wall. We have taken the same weight over 
the linear sum since there is nothing to prefer one m from 
another. Therefore, the behavior of J-i^(a~^z) would be 
of the form 

Ci(m-Vff)i/2+«^^ + C2(m-VH)i/2-^^ , (25) 

where Ci , C2 are some constants which are not important 
for our present purpose. The boundary condition (|24p 
then implies 

C(l/2 + i£;)Cir^/'+'^ + C(l/2 - iE)C2r^jl^'"^ = 0(26) 

The scalar field eigenvalues which are independent of 
short distance cutoff can be obtained from (I^S)) by sep- 
arately making the coefficients zero, which implies that 
eigenvalue equation for the scalar field is 

C(l/2 + z£;) = 0. (27) 

The scalar field eigenvalues thus coincide with the non- 
trivial zeros of the Riemann zeta function. 

V. Discussion. — Although we considered the back- 
ground to be a (3 -|- l)-dimensional Schwarzschild space- 
time, the analysis is applicable to a more generic black 
hole spacetime, because the effective Hamiltonian for the 
field in the near-horizon region for a class of black holes 
are known to have the scale invariant form (Tsj . 

Note that in a similar discrete scaling was made 
and a sum taken over the scaled states of the sym- 
metrized operator D = xp. Then the zeta function 
C(l/2 -I- iE) came out as the overall coefficient of the 
resulting state. Unlike in our case, at the zeros of the 
zeta function, therefore, the wave function ^(x) was a 
null state for all x. 
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